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ABSTRACT 


This report describes ongoing activities performed at the Vir- 
ginia Institute of Marine Science (VIMS) for NASA Wallops Station to 
study the applications of remote sensing in marine science. During 
1973 scientists at VIMS utilized remote sensing in three programs: 
tonal variations in imagery of wetlands, use of the thermal infrared 
to delineate the discharge cooling water at the Virginia Electric and 
Power Company (VEPCO.) nuclear power station on the James River, and 
the use of aerial photography to determine the volume storage function 
for water in the ■ marsh-bay complex fed by Wachapreague Inlet on the 
Eastern Shore of Virginia. 

One program was to determine the causes of tonal variations -in 
wetlands imagery and delimit their scope. Results from a predominately 
single -species salt marsh show variations in tone large enough to 
cause confusion in an automatic mapping scheme or bioproductivity 
studies. Film density-bioproductivity correlations for a given sample 
exhibit variations at different points on the film frame. 

Study of the cooling water discharge from the VEPCO nuclear 
power plant and its ef fect on the local marine environment is 'part of 
a continuing program funded by NASA Wallops and the Atomic Energy Com- 
mission. Thermal infrared imagery from a scanner operated by NASA 
Wallops was used to obtain an instantaneous surface thermal map of the 
area. Rectification of the strip film from the scanner, thermal- 
spatial resolution, and the correlation of the geometrically corrected 
scanner output with ground recorded information is discussed. Some 
insights into the plume behavior are presented. 
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Finally, a new application of remote sensing was successfully 
tested using aerial photography to determine the tidal volumes for 
various tidal stages within a marsh-lagoon system in a harrier island 
complex.. The capability of black and white .infrared film to distinguish 
a water-land interface was used as a contouring machine to determine 
the area within the system covered with water, and the volume storage 
function . 


iv 



INTRODUCTION 


For the past three years the Virginia Institute of Marine Science 
and NASA Wallops Station have worked together in a program to apply 
remote sensing in marine science. VIMS has studied the southern Chesa- 
peake Bay and, its tributaries extensively over the past twenty years, 
with most of the data obtained in the field using direct sampling methods. 
VIMS and NASA have recognized the fact that some of the necessary in- 
formation can be obtained more easily using remote sensing (and exist- 
ing fields of study thereby complemented). Several areas of marine 
research have consequently been investigated for the application of 
remote sensing. This and two earlier annual contract reports (1971, 

1972) detail the progress ache ived. 

One of the initial programs utilizing remote sensing was the 
remote detection of oil on water, and a study of oil slick behavior. A 
great deal was learned using different film and filter combinations as 
well as the thermal infrared. With most of the experimental releases 
made in the ocean we believe this work summarized well the best choices 
available using conventional aerial photography, and the advantages 
and limitation of the thermal infrared. 

A second area where remote sensing has proven valuable is a 
study of barrier island inlet dynamics conducted on Virginia T s Eastern 
Shore at Wachapreague Inlet. Here remote sensing complemented extensive 
bathymetry and current recording by giving a plan view of changing sur- 
face features. 

A new application of remote sensing to barrier island inlet 
dynamics is presented here. A series of repetitive overflights using 



panchromatic infrared film was executed to quantify the water volume 
storage function of the marsh-bay system fed by .the inlet. The method 
and results of the work are discussed in this report. 

The study of tonal variations in wetlands imagery is based upon 
a phenomenon recognized and studied in agricultural crops. Tonal vari- 
ations have been shown to cause significant confusion in the inter- 
pretation of wetland photography (see Gordon and Penney, 1973), but 
have not been studied in detail. Analysis this year has focused on 
determining- the magnitude of the variations, and possible correlation 
with measured plant parameters for a homogeneous salt marsh. 

Use of thermal infrared imagery at Hog Point is the continuation 
of an ongoing program of extensive ground measurements and aerial 
thermal mapping to delineate volume and surface thermal plume from 
the Virginia Electric and Power Company (VEPCO) nuclear power plant. 

Scientists at VIMS who have utilized the remote sensing avail- 
able from NASA Wallops have been pleased with the information they 
have obtained and the extra dimension it has added to ongoing programs. 
During the past three years of this contract there have been important 
uses of remote sensing as well as new techniques for obtaining support 
data. We hope those who read this report and the past two yearly re- 
ports will find information of benefit to their ongoing or projected 
research . 
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TONAL VARIATIONS IN WETLANDS IMAGERY 


Introduction 

Even a nonexhaustive review of wetlands research literature 
reveals at least three types of research problems of management 
Interest: (1) defining wet land- dryland boundaries; (2) mapping and 

determining areas of species or communities and (3) determining species 
or community productivities. Egan and Hair (1971), Wobber and Anderson 
(1972), Gallagher, Reimold and Thompson (1972), and Klemas, Deiber, 
Bartlett, Crichton and Fornes (1973), have approached the first two 
problems using various remote sensing media. Analytical techniques 
have ranged from subjective comparisons of color and tone to auto- 
matic densitometry and planimetry. The time consuming nature of 
interpretive studies is universally recognized and some mention is 
often made of tonal or reflectance variations as presenting inter- 
pretive difficulties. For instance, Egan and Hair determined delimit- 
ing tonal ranges for wet and drylands over the film frame for their 
geographic area. Klemas, et al. point out that n The uniform flatness 
of marsh topography eliminates variations in reflectance due to slop- 
ing surfaces and shadows.” Gallagher, et al. recognize, to a degree, 
that viewing angle changes are akin to topographical angle changes and 
use the variable illumination in a Kelsh type plotter to balance out 
tonal variations in conjugate areas of stereo pairs. 

Productivity estimates using correlations between measured film 
densities or reflectances and plant parameters have been mentioned by 
Gallagher, et al . To date, they have presented only subjective corre- 
lations and no detailed quantitative density measurements. Their 
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feeling that quantitative measurements must always be complemented by 
some degree of "marsh awareness” or interpretive ability is in some 
degree supported by this work. Previous work at this institution, 
Gordon and Penney (19 72), has demonstrated that when fT marsh awareness” 
is purposely ignored, tonal variations caused up to forty percent 
changes in area changes of "plant” or scene types for conjugate areas 
of the film. 

Hoffer, Holmes and Shay (1966), Egbert and Ulaby (1972), Steiner 
and Haefner (1965), Suits (1972) and Silvestro (1968) have made quanti- 
tative and theoretical investigations into the relationships between 
■scene type, tone, reflectance, and viewing geometry. In general, image 
tone is a complicated, but determinable, function of incident solar 
flux, angular dependent reflected energy, atmospheric effects, camera 
radial fall-off function, and film (type and processing effects). 

Concerning ourselves only with one mission over a short time, 
tonal variations are a less complicated function of scene type, 
radial fall-off, and changes in scene reflectance with viewing geometry 
For one scene type, in terms of light reflecting parameters, variations 
in reflectance or density can be easily determined using procedures 
detailed by Silvestro (1968). Using these techniques quantitative 
density or reflectance can be correlated with plant parameters com- 
monly measured by ecologists. 

Technique 
Field Work 

Field work was conducted in a high salinity Spartina alterna- 
flora marsh known to contain three growth forms of Spartina and small 
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amounts of a Salicornia species. Quarter meter square clip plots 
were gathered in thirty-three locations predetermined using aerial 
photographs. An attempt was made to sample all scene types adequately 
and to place the images at many locations on the composite film frame. 
The grass samples were analyzed for stem density, length, and diameter 
then dried and weighed. A summary of measurements for 13 samples used 
for analysis is presented in Table I. 

Histograms for height, weight, stem density, and basal (or 
total cross-sectional) area revealed that, the traditionally recog- 
nised growth forms — tall, short and intermediate — are not completely 
distinguished by these parameters. For instance, "tail” and ,f short” 
forms had similar weight and basal area ranges, but different height 
and stem density ranges. The "intermediate” form intergraded for all 
parameters . 

Remote Sensing Photography 

On November 20, 1972, a NASA-Wallops Station remote sensing 
aircraft (C-54), configured with a T-ll mapping camera (black and 
white infrared film - 89B filter) and an I 2 S four band camera made 
four passes over the marsh at -1200 EST (low tide) and four at 
— 1400 EST (low tide +2 hours). All runs were at 5,000 feet (scale 
1,10,000) altitude. Two additional low altitude passes were made 
over a set of five different gray panels located adjacent to the 
marsh and used to facilitate reflectance calculations. 
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TABLE I 


Sample # 

AV. Dia . 
(mm) 

Av . Stem 
Length 

Total 
Dry Wgt . 
(grams) 

# stems 

Approx o 
Basal area 
(mm 2 /plot) 

- 15 

3.5 

30.8 

38.5 

134 

1146 

33 

3.0 

38.7 

31.3 

72 

509 

18 

3.8 

33.4 

25.1 

49 

556 

19 

3.6 

43.8 

75.0 

125 

1272 

25 

2.4 

29.6 

31.4 

28 

127 

11 

6.4 

90.3 

101.4 

39 

1255 

9 

■ 9.9 

115.1 

182.1 

22 

1693 

14 

5.2 

60.5 

47.1 

30 

1146 

24 

4.2 

47.6 

100.2 

133 

1843 

10 

2.9 1 

23.9 

44.3 

108 - 

2739 


3.2 s 

20 

56.3 

252 

(combined) 

22 

5.8 

82.6 

115.7 

54 

1427 

26 

3.0 

30.9 

71.0 

201 

1421 

28 

3.0 1 

. 33.6 

104.6 

204 

2993 


5o0 2 

19 

12.8 

79 

(combined) 


1 Spartina 
3 Salicornia 
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Calibration 


In order to separate density changes caused by the sensor 
system from those related to the scene and intervening atmosphere it 
was necessary to calibrate the cameras. To correct the mission film 
for non-uniform camera- transmission effects (which we summarize 
"radial fall-off") a special calibration source was employed (jointly 
designed and constructed by Dr. J.D. Oberholtzer and Mr. Charles 
Vaughn of NASA-Wallops and the authors). The calibration source is, 
ideally, placed in contact with the camera lens cone (Figure 1) and 
the camera cycled in the normal manner. Since the source has a uni- 
form energy distribution over an FOV greater than the camera, off-axis 
tonal changes can be determined. These were found to be radially 
symmetric for the T-ll (Figure 2) and asymmetric for the I 2 S (Figure 3) 
For the T-ll the' relative exposure between the center and an off-axis 
point is E(0)/E(O), where 0 is polar view angle (PVA)(see Figure 4). 
Since film density, D, is a logarithmic function of exposure, E, the 
normalized (to the center) density correction is AD = D(0) - D(o). 
Densities at any point on a mission frame can then be corrected when 
the PVA of the point is determined. In this experiment, the radial 
fall-off calibrations were exposed and processed separately from 
actual mission film. Both pieces of film were normalized to a gamma 
of one. The asymmetry of the I 3 s camera system is unfortunate and 
may have been caused by extraneous light leaks due to lack of contact 
between lens cone and calibration source, or flare, which the I 3 S 
occasionally exhibits on normal missions. The asymmetry precluded 
easy correction of’ radial fall-off. 
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Figure 1. Diffuse source for camera calibration. 
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Figure 2. Density correction factors for diffuse calibration 
of radial fall-off function. 


CO 

VC 


T-11 


I. 


(T 



10 20 


— 1 1 

40 50 


EW ANGLE (PVA) 


9 . 


i 

1 I 

i 

j 

i \ 

i 

Figure 3. Density values- from an I 2 S~ negative exposed with j. 

the diffuse source . 




Figure 4. Viewing geometry 
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Results 


^positional Analysis 

Figure 5 is a smoothed histogram of fall-off corrected mean 
densities, from a positive, for the three traditional plant types and 
mud flat from the 1200 EST' flight, IR band . Obviously, the three types 
would be frequently mistaken for each other when viewing geometry and 
tT marsh awareness” are ignored. Plots of mean corrected film density 
versus stem density (Figure 6), height (Figure 7), weight (Figure 8), 
basal area at 1200 EST (Figure 9), and 1400 EST (Figure 10) show no 
obvious functional relationships (numbers refer to Table I samples). 

For example, images (of plant samples) having a contrast of approxi- 
mately 0.2 density units at all (present) film positions were found 
to occupy overlapping stem density, height, and weight ranges. So, 
if viewing geometry had been ignored, a contradicting, false inter- 
pretation would have been made about film density parameter corre- 
lations. Basal area, which is inversely proportional to the soil 
visible, gives much less contradictory results, even ignoring viewing 
geometry. Above a basal area of approximately 1400 mm 2 /plot there is 
a sharp drop in mean film density (or increase in reflectance — see ap- 
pendix). This density or reflectance change is significant at the 
10% level (two sample t test) if the groups greater and less than 
1400 mm 8 /plot are compared. The density difference is less startling, 
though still significant at the 10% level, on the 1400 EST flight for 
reasons not presently clear to us, but probably involving a change in 
the geometry dependent reflectance function at the lower sun angle. Com- 
plete analysis of such an effect requires a more extensive sampling pro- 
gram with different basal area groups moving through the same point loci. 
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PER CENT SAMPLE 


Figure 5* Film density histogram. 
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Figure 6. Corrected mean film density vs. 
T-ll IR band (~12Q0EST). 



PLANT DENSITY (number of 


stem density 


28 |p 

— I 1 *“T 1 

280 320 360 400 

stems /plot) 


20 40 60 80 100 120 
AVERAGE HEIGHT (cm) 




Figure 8 . 
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Figure 9.* Mean film density vs. basal area plot T-ll IR band 
(-1200 EST). 
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Figure 10. Corrected mean film density vs. 
IR band (1400 EST). 
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Positional Analysis 

Figures 11, 12, 13 and 14 are various ways of presenting den- 
sity data positionally- The isoline plots in 11 and 12, contoured by 
someone unfamiliar with tonal variation concepts, are for the < 1400 mm / 
plot group. The patterns, are similar for the two flights, having 
slightly greater rates. of density change on the northerly side of the 
frame (or sun azimuth line) than the southerly side. East-west dif- 
ferences may be indicative of some real functions of viewing geometry, or 
simply, sampling inadequacies. Further work in this type of marsh 
could lead- to generation of a reliable density correction mask for view 
angle effects. Ideally, a mask for each of the important marsh species 
would make productivity mapping more accurate and help identify view- 
ing geometry regions in which species discrimination would be more 
accurate . 

Figure 13 presents the <1400 mm 2 /plot with the composite film 
frame divided into, columns perpendicular to the flight line (and point 
loci). This approach divides the total sample . into two azimuthal groups. 
Points between 0° to 90° AVA and 2 70° to 0° AVA are one group,'. Points 

* . i . > % 

90° to 270° are the other. Column -4 to +4 represent equidistant di- 
visions of a sample point locus from its first entry into the frame 
to its exit point. Column 0 is the closest approach to the frame 
center. It is interesting that between group contrast changes slightly 
over the frame, but, as in the isoline presentation, rate of density 
change appears to be greatest on the north side of the frame. The area 
in black is the overlap of the one standard deviation units for each 
group by column. Variability in this overlap may be due to sampling 
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Figure 11. Contoured composite photo for samples <1400 mm 3 / 
plot (~1200 EST). 
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Figure 12. Contoured composite photo for samples < 1400 mm 3 / 
plot (~1400 EST). 




Figure 13. Composite frame points grouped by columns perpen- 
dicular to flight line . 
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deficiencies. Figure 14 is the composite frame divided into polar 
view angle annuli. Differences in contrast and standard deviation 
overlap may again be sampling deficiencies. Essentially, the column 
approach ignores polar view angle and the polar view angle approach 
ignores azimuthal view angle. 

An analysis of variance would require different sampling, 
but it seems clear that "frontlit" samples have higher reflectance 
than "backlit", and that reflectance increases (density decreases 
on a positive) with increasing polar view angle. Both conclusions 
are supported by work in agricultural remote sensing literature. 

Raw-Corrected Data 

Figure 9, from T-ll data, basal area at 1200 EST, compares 
fall-off corrected and uncorrected data. It is important to note 
that the reflectance increase at 1400 mm s /plot would not have been 
visually noticed in the uncorrected data. plot. In fact, a two sample 
t test on the raw data for greater than 1400 mm 3 /plot and less than 
1400 mm s /plot sample groupings is not significant at the 10 % level. 

The mean T-ll corrections were used on the I 2 S red and green band 
data, making the very risky assumption that the I 2 S fall-off is indeed 
symmetrical and of the same general shape as the T-ll. Figures 15 
and 16 show red and green band raw and corrected densities plotted 
with basal area. Possible fitted lines are shown with dashes. The 
data is suggestive of no correlation In the red and slight correlation 
in the green- These differences between raw and corrected data imply 
that caution must be used when interpreting quantitative, but un- 
corrected, data. 
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MEAN FILM DENSITY 


Figure 14. Composite frame points grouped in polar view angle 
annuli . 
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Figure 15. Mean film density - raw & "corrected” - basal area 
I s S red band (~1200 EST>. 
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figure 16. Mean film density,- raw 6 "corrected” - basal area 
I 2 S green band 1200 EST) . ^ 
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Reflectances 


It is possible, following Silvestro T s method (1968) to calcu- 
late reflectance of scenes from corrected densities. Film exposure 
on the negative is proportional to reflectance in an equation of the 
form 

E (X) = oc (X) r (X) + B (X) 

where a(X) is a constant (over the exposure time but dependent on 
solar flux) and B(\) is that energy from atmospheric scattering 
added to the imaging light. In our case X represents the 0.7 - 0.9 
micrometer band (B &'W IR) . Three such equations of this sort were 
written for two ground panels of known reflectance as well as for the 
scene. Solving the three equations for target reflectance, r (X), the 
following is obtained: 

„ _ [ExU) - E B a)] r w CX) + [E w a> - E x (X>] r R CX) 

x E W (X) - E B (xr 

where E w , r , Eg, r^ are exposure and reflectance for the white and 
black ground panels, respectively. 

Reflectance data for the IR band is presented in the appendix. 
Plotted with plant parameters, reflectance showed essentially the 
same behavior as mean corrected densities. The 1400 mm 2 /plot dis- 
continuity was still evident. 

Conclusion 

It seems clear why, when marsh awareness is ignored, subjective 
color or tone keys can prove inconsistent or ineffective. Film den- 
sity in all bands varies due to lens illuminance (radial) fall-off, 
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and scene viewing geometry. On a positive, with increasing polar 
viewing angle, marsh scene reflectance increases and therefore den- 
sity decreases. Lens illuminance fall-off has the opposite, in a 
sense,, partially balancing -effect. As a result of these opposite 
trends in density change, use of an anti- vignetting filter causes 
tonal variations due to salt marsh scene viewing geometry to be 
more visually noticeable. Analysis of some six hundred density measure 
ments presented here showed that TT raw Tf illuminance corrected measure- 
ments can give different results when plotted with various plant 
parameters. Plant samples at the extreme ends of the height, stand- 
ing crop biomass, and stem density ranges often had closely similar 
mean (corrected or raw) film densities. For this data set, there 
was a statistically significant change in mean corrected film density 
at a basal area of ~ 1400 mm 2 / quarter meter 2 sample plot. Further 
work, in progress, may or may not make this conclusion generally ap- 
plicable for this type of marsh community. However, the implication 
is that in the near IR band reflectance is not as dependent on 
biomass, stem density, or height as it is on amount of soil visible. 

Reflectance (or density) plant parameter correlations will con- 
tinue to be investigated^ by the authors. An attempt is being made to 
compare only samples having the same viewing geometry. Use of color 
IR film will speed multiband data analysis by eliminating the need 
to locate each image point more than once in the densitometer.. Ground 
panels and densitometry will continue to be used, and fall-off cali- 
brations will be performed on the mission film to eliminate gamma 
normalization between films. 
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HOG POINT THERMAL STUDIES 

Introduct ion 

Thermal studies of the James River in the vicinity of the Vir- 
ginia Electric and Power Company (VEPCO) nuclear power generating 
plant (Surry County, Virginia) have continued during 1973. The work 
done by VIMS under contract to the Atomic Energy Commission (AEC) and 
NASA Wallops this past year has been directed toward monitoring the 
immediate section of the river near Hog Point throughout the initial 
period of operation of the power plant. The background conditions of 
the environment have been studied during 1971 and 1972 with an in- 
strumented boat, in situ tower, and airborne thermal scanner (see Bolus, 
et al., 1971; Shearls, et al., 1973; Gordon and Penney, 1972)-. 

Research under the NASA Wallops contract has been directed 
toward the use of the thermal infrared, its interpretation, and appli- 
cation to the problem of mapping the temperature structure in a time 
dependant estuary. Described in this section of the report are the 
scheme used for rectification of the infrared imagery, determining 
the resolution of the scanner, and results of a thermal flight during 
high slack tide and the beginning of ebb flow in the outfall vicinity 
of the James River with the nuclear power plant in operation. 

Thermal- Spatial Resolution 

In order to locate discrete positions on the infrared scanner 
film record such as the instrumented boat or the data acquisition 
tower it is necessary to have an object which emits either more or 
less energy than the surrounding water and is within the spatial. 
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resolution capability of the scanner. From past experience it was 
found that the 30 foot instrumented boat could not be located on the. 
film at the working altitude of 2,000 feet. An experiment was devised 
to run. concurrently with an actual data flight in the field. To get 
an order of magnitude measurement .of the thermal- spatial resolution, 
overflights were requested at 500, 1,000, 2,000 and 5,000 feet alti- 
tudes over a boat anchored at the tower with a special infrared target. 
To make the target two sheets of masonite x 4 f x 4 T ) were painted 
silver on one side (low emi.ssivity) and flat black on the other side 
(high emissivity) . The energy source was the sun and the panels were 
oriented during the day to form a 4 f x & 8 T blackbody (both large 
and small emissivity). A small hole was drilled in the edge at each 
panel to accept a special thermistor for measuring the internal temp- 
erature. Figure 17 is an unsmoothed plot of the temperature variation, 
for each configuration. The fluctuation in values, especially noted in 
the black panels, was due to interruption of sunlight by clouds. Be- 
tween 1500 and 1600 hours the sun was completely obscured and the black 
panel reached an equilibrium temperature of 20-21° C from indirect 
radiation. Results of the experiment are summarized below in Table II. 

Table II 


Altitude Black 
4x8 


Silver Black 

4x8 4x4 


.5 X V V 

IX V M 

2 X M N 

5 X N N 


V 

M 

N 

N 


V - visible, M- marginal, N - not visible 
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To be classed as visible an observer must be able to recognize 
the image on the scanner film placed on" a light table with no magni- 
fication. Use of this type of marker system is, therefore, not 
feasible for work done at altitudes greater than 1,000 feet. A sug- 
gestion we would like to try during the course of next year T s experi- 
ments is use of a charcoal brazier in a floating platform (H. Jones, 
Photo Science, Inc., personal communication). Although spatially 
smaller than one of our panels, the. charcoal at its very high temp- 
erature should radiate enough energy to image as a discrete spot on 
the scanner film. 

4 

Geometric Rectification of IR Scanner Film 

The utility of strip film from an infrared scanner is greatly 
reduced if features which are thermally imaged cannot be put into 
proper perspective and correlated with other sources of information. 

The problem of geometric distortion is particularly acute over water 
where there are relatively few stationary features to help orient the 
thermal image and remove the inherent distortion. Causes and mani- 
festations of distortion in IR strip imagery are discussed in articles 
by Sabins (1973) and Masry and Gibbons (1973). In order to use the 
imagery to outline the thermal plume in its proper perspective the 
progressive compression of information toward the film edges caused 
by constant angular scanning and linear recording must be eliminated. 
One method of accomplishing this is to use the film r s center portion of 
roughly 40° coverage where distortion is minimal, and by flying para- 
llel lines to gain the desired lateral coverage. Reconstructing the 
overlapping strips becomes very difficult over water, however, for 
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reasons stated above. Another solution is to reduce the scale of the 
imagery, thereby confining the area of interest to the central portion 
of the strip. This is possible only to the point where too much in- 
formation from the scene is lost through the atmosphere. As a result 
of experimentation, 2,000 feet was chosen as an optimum altitude for 
recording the maximum amount of information at the smallest scale 
over the Hog Point site. 

At an altitude of 2,000 feet the thermal plume is generally 
confined to one strip width flown such that the shoreline is at the 
very edge of the film format. In this situation rectification is 
needed to accurately locate the plume and determine its true shape. 

In order to accomplish this a simple mapping transformation was 
employed. A rectangular grid was used to overlay the scanner film 
and a special grid with progressively wider spaced lines in the cross- 
track dimension (symmetrical about the nadir) was drawn. A point by 
point transfer of features of interest from the film to the special 
grid served to "undistort” the scanner imagery. The method was 
tested using imagery flown east- west over Hog Point at the power 
plant site and linear features such as roads and power line clearings 
were checked against a recent map of the area made from aerial 
photography. Correlation was in general good although there Was 
some misplacement of features caused by aircraft crab, altitude change 
over the run, and speed variation. Still, agreement with the .base 
map was better than 10% (displacement of known features from their 
true location) . 
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An improvement of the rectification scheme, to be tried next 
year, will rely on automation. Imagery will be digitized on the NASA 
Wallops scanning microdensitometer, the rectification performed on 
the magnetic tape with computer processing, and the resulting information 
used on the same densitometer to expose a new piece of film. The new 
film will contain the original radiometric information relatively free 
of geometric distortion. 

Thermal Flight 4/17/73 

On April 17, 1973 the NASA Wallops C54 equipped with the AAD-2 
infrared scanner and a T-ll mapping camera (S0397 color film) made 
repetitive flights over the VEPCO nuclear power plant discharge jetty 
and adjoining section of the James River. Ground truth during this 
period was gathered from the instrumented boat running transects, a 
second boat anchored at tower #6, and some instrumentation on. the 
tower itself. Measurements made by the moving boat were water temp- 
erature at 0.5, 3 and 6 feet below the surface, air temperature at 
3 and 6 feet above the surface, dew point temperature, and salinity. 

The second boat, in radio communication with the aircraft, recorded 
water -temperature during the day at the 0.5 foot depth (as well as 
performing the spatial-temperature resolution experiment). Since the 
Braincon data acquisition system was not installed on the tower, . wind 
velocity only was recorded. 

The aircraft made four groups of passes over the area beginning 
at 1040, 1200, 1330, 1530 EST • Each group consisted of three passes 
perpendicular to the discharge shoreline (center line over the power 
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plant plus one line on either side) at 500 feet altitude, and four 
passes parallel to the shoreline^ and directly over tower #6 arid the 
stationary boat at altitudes of 500, 1,000, 2,000, 5,000 feet. The 
instrumented boat began its run on transect #4 (T4) at 1102 EST and 
completed a partial run through transect #8 (T8j at 1211 EST* Engine 
trouble precluded making an entire run or repeating the transects as 
had originally been planned. The stationary boat was put on location 
at the tower in the morning (1000 EST) and remained until 1600 when 
the aircraft made its last pass. Before returning to the dock the 
boat took a final temperature reading approximately 100 feet outside 
the discharge jetty. 

As a background for the following discussion Table III is in- 
cluded to summarize the status of factors important to the plume behavior 
Figures 18, 19, 20, and 21 show the rectified thermal plume outline 
on a map of Hog Point. All four are to the same scale and each made 
from one strip of scanner imagery flown from north to south at 2,000 
feet altitude parallel to shore and directly over tower #6. Figure 22 
is a composite made from four of the original 70 mm scanner strips 
taken at 1140, 1240, 1440, and 1600 EST. The two dark lines running 
longitudinally either side of the center are artifacts caused by the 
recording process and should not be interpreted as a part of the scene. 
Figures 18 and 19 show the transects followed by the instrumented boat. 
The temperature contours are drawn from the data taken by the boat at 
the 0.5 foot depth. Refer to Appendix B for the original data re- 
corded by the boat on April 17. Figures 20 and 21 have no boat transects 
or temperature contours since the overflights were made after the boat 
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INAUTICAL MILE 

Figure 18. Rectified map showing scanner derived thermal surface 
plume (1140 EST), boat transects with times (EST), 
and boat derived isotherms at 0.5 foot depth. 



hog point 



I NAUTICAL MILE 

Figure 19. Rectified map showing scanner derived thermal surface 
plume (1240 EST), boat transects with times (EST), 
and boat derived isotherms at 0,5 foot depth. 
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Figure 21. Rectified map' showing scanner derived thermal surface 
plume (1600 EST). 




Figure 22. 


Composite made from four original 70 mm infrared scanner strips 
(8-14 (jm). Times from top to bottom - 1140, 1240, 1410, 1600 EST • 
Warmer objects appear lighter. Altitude was 2,000 feet over 
outfall jetty. 





runs were finished. Figure 23 is an unsmoothed plot of temperature 
at the tower throughout the day at the 0.5 foot depth. Finally Fig- 
ure 24 is presented to show the discrete areas of differing surface 
temperature as recorded by the IR scanner. Shown are four black and 
white reproductions of four color encoded, density sliced portrayals 
of the scanner strips in Figure 22 (note only a small portion was used). 
Temperature of the surface layer of water increases from region A to 
B to C. Variations or banding within each region is caused by shades 
of the original three colors (blue, violet, brown) which further sub- 
divide the major regions into subregions* For purposes of discussion 
in this report only the major regions (A, B, C) will be used. The 
shoreline and discharge jetty have been encoded black for orientation only 


Table III 

Conditions at Hog Point on April 17, 1974 
Tide (predicted) 

Max High Max Low 

Slack Flood Water Slack Ebb Water 

0803 1017 1111 1301 1634 1726 

Wind 

0-5 knots SW at 1200 

Nuclear Plant Operational Parameters 1 

Intake Temperature 3 57.8°F(14. 3°C) 

Discharge Temperature 66 .4°F(19 .1°C) 

Operation Unit 1 Unit 2 

86. 8% 3 43.2% 

1 From VEPCO records 

2 At 1200 EST 

3 Based on a 100% output of 788 megawatts and 11.9 x 10 9 BTU/hr 
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Figure 24. Black and white print of color encoded representations of 

original scanner imagery shown in Figure 22 (only a portion 
was used, and scale is not the same as Figure 22, or original). 
Radiometric temperature increases from blue to violet to 
brown and from darker to lighter shades (A to B to C) . 
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By analyzing these figures and material in Appendix B the fol- 
lowing conclusions were drawn- The temperature at the tower as recorded 
by the stationary boat and the temperature recorded by the moving boat 
at the time of its closest approach to the tower agree to within 0-l°C 
The overall shape of the plume as recorded by the thermal scanner 
agrees very well with the tide stage in Table III. Figure 22 shows 
the surface plume -at 1240, very near the time of predicted slack water, 
to be flowing straight from the exit canal jetty with slight curvature 
in the direction of ebb flow (down river toward Hog Point). In the 
morning (1140) the surface plume was moving into Cobham Bay and forming 
a well defined envelope extending approximately 1.5 miles to the north 
of the outfall. The presence of isolated warm spots within the plume 
is evident in Figure 24. Figure 18 (the rectified map with boat tran- 
sects taken at 1140) shows a degree of correlation between the inner 
warm plume and the 14.0°C contour at the 0.5 foot depth. Care should be 
exercised in interpreting the rectified maps with boat derived temp- 
erature contours since these are drawn using data taken at different 
locations separated by up to one hour in time. By the 1408 pass, the 
rectified map (Figure 20), the scanner reporduction (Figure 22), and 
the color encoded representation (Figure 24) all show a warm "finger” of 
surface water beginning to move toward the point (downriver), surrounded 
by a new well-defined surface plume envelope approximately 0.5 by 1 
mile. Remnants of the earlier envelope can still be seen in Cobham Bay 
in Figures 22 and 24. There is good correlation between the warm 
part of the plume (which is very close to the tower location) and the 
2°C rise in water temperature recorded by the stationary boat between 
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1330 and 1400. When the 1410 pass was made the warm "finger” had 
apparently . just left the tower as evidenced by its location in Figure 20 
and the corresponding drop in water temperature at the 0.5 foot level. 

By 1600, near the time of maximum ebb flow, the scanner shows the 
plume to be moving from the jetty, folding close to shore, and merg- 
ing with other warmer water from upriver rounding Hog Point (refer to 
Figures 21, 22, and 24). The sharp bend made by the plume leaving 
the jetty indicates a water current which would preclude the formation 
of a well-defined surface plume as noted earlier. 

An attempt was made to correlate the color encoded surface plumes 
in Figure 24 with temperature measured at the 0.5 foot depth. To test 
the validity of the correlation we hypothesise it is possible and ex- 
amine our data for supporting facts. The following assumptions are 
made: the density of the IR film image varies linearly with surface 

temperature for equal emissivity targets (over a small temperature range), 
the color encoding process assigns colors and shades of colors- to equal 
density ranges, and the increase in path length for images at the edge 
of the scanner field of view does not cause a significant energy loss 
at 2,000 feet altitude. By the hypothesis the brown shades (c) in 
Figure 24 will be assigned the range of their highest and lowest measured 
temperatures, 16.5°C to 19.0°C. 16.5°C was measured at the tower at 

1400 when the warmest part of the plume surrounded the tower; 18.6° was 
measured by the boat at 1615 in the immediate vicinity of the jetty while 
VEPCO lists a temperature of 19.1°C at noon within the jetty. Using the 
assumptions, all shades of violet (B) will represent a 2.5°C range as 
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will all shades of blue (A)- The temperature ranges assigned to the 
colors would then be 11.5 - 13.9°C for blue, 14.0 - 16. 4° C for violet, 

and 16.5 - 19°C for brown. Referring to Figure 18 and the color en- 

/ 

coded representation (Figure 24) the instrumented boat passed through 
the area of the plume coded violet at the -time the scanner was making 
the 1140 overpass. Temperatures recorded by the boat at 0.5 foot 
depth are in the range of 13-14° C which would be coded blue using the 
hypothesis. Shortly thereafter the boat passed out of the plume, 
coded violet, and into the area coded blue in Figure 24. However, 
toward the end of transect #7 and beginning of #8 there are temperatures 
at the 0.5 foot depth in the 14-15°C range which would be coded violet 
using the hypothesis. 

The total temperature range at the tower was 2.5°C while that 
recorded by the moving boat was 3°C. Taken together this represents 
a range of 4.5°C, whereas the scanner responded to a hypothetical 
7.5°C temperature spread. The above facts indicate that the surface 
thermal plume will not always have a temperature the same as that 0.5 
feet below the surface. Indeed, boat data from depths of 0.5 and 3 
feet indicate that the surface layer of water may be several degrees 
warmer than the subsurface in some locations and the same in other 
places. The original hypothesis is not valid. We are able to define 
the extent of the surface plume and to delineate the warmer and cooler 
areas, as well as the maximum plume temperature of 18-19°C. Temperature 
contours derived from the infrared scanner data are not possible, how- 
ever, without additional ground truth or scanner calibration. 
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Continuing Studies 


The concentration of work for 1974 will be on correlating in- 
frared imagery and ground data with the power plant in operation. 
Overflights will be scheduled during periods of extreme and typical 
plume tidal excursion, and during times of the year when fresh water 
discharge is maximum or minimum. We suggest, from past experience, 
planning extra overflights in order to account for weather problems, 
equipment failure, or other events which might, arise. One flight 
per month, beginning in March when the instrumented boat again resumes 
its schedule of bi-weekly data runs, should be sufficient. 

Coincident with the thermal flights the experiment to find a 
suitable target easily recognized on the thermal imagery will con- 
tunue using a pan of hot charcoal anchored in a special floating plat- 
form. It is hoped the temperature of the charcoal will offset its 
limited area enough to produce a white spot on the film positive. 

One additional flightline, parallel to the existing north- south 
line and slightly farther offshore, should be added to image the 
tower toward the frame edge to determine how much energy is lost due to 
the increase in atmospheric path length. This would also allow more 
reliable placement of the outer edge of the plume and give a better 
level for the background surface water temperature near the river 
channel . 

A major part of the effort this year will be to examine the 
surface temperature of the water in relation to the temperature at 
the 0.5 foot depth to determine the extent to which the plume,’ as 
imaged by the scanner, is a surface layer of warmer water over cooler 
subsurface water. Use of the floating thermistor employed in the 
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oil experiments conducted by VIMS and NASA Wallops (Gordon and Penney 
1971 and 1972) as well as a radiation thermometer in the aircraft to 
better assign temperature to scanner film density is suggested. 

Finally, a further improvement over last year f s techniques 
will be automation of the geometric rectification-density enhance- 
ment process using the NASA Wallops scanning microdensitometer and 
computer facilities as described previously in this report. 



DETERMINATION OF TIDAL VOLUMES IN BARRIER ISLAND SYSTEMS 
BY SEQUENTIAL BLACK AND WHITE IR IMAGERY 

I Introduction 

Much of the Atlantic and Gulf coasts are characterized by 
barrier island complexes comprised of marsh- lagoon systems subject 
to tidal flushing by inlets. These coastal areas are experiencing 
increasing pressures of development by man and as a result the water 
quality of these ecosystems is frequently degraded. Enlightened 
management of these resources demands a rather complete understanding 
of the flushing characteristics of the basin and interchange of water 
via the inlet. The stability of the inlets themselves Is central to 
flushing of the interior basins since empirical evidence (O'Brien, 
1969) indicates that the cross-sectional area of the inlet channel 
controls the tidal prism gated by the inlet. The inlet channel area, 
moreover, is controlled by the magnitude of littoral drift (tending 
to choke the inlet) and the hydraulic currents generated by the 
difference of water levels due to tide fluctuations tendingr to open 
the inlet. Thus the inlets can be viewed as the keystone of the 
flushing of the barrier island system as well as access points to the 
interior by boats. 

The dynamics of inlets are very incompletely understood. It 
Is not possible, for example, to predict specifically how an 'inlet 
will change under given conditions of littoral drift tending to close 
the channel. Recent research on inlet dynamics, however, tends to 
confirm the fact that the volume of water gated by the inlet to the 
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interior controls the overall stability of the channel. Furthermore 
the storage characteristics of the basins (volume stored as a function 
of water level) determines the distortion of the dis charge -time relation 
ship at the inlet channel which thereby determines whether or, not sand 
will be advected into the interior of the inlet channel (Mota Oliveira, 
19 70 ). 

The Problem . The determination of tidal volumes for various 
tidal stages in the basins behind barrier island systems is essential 
to : 

1) understanding the behavior of and the design of tidal 
inlets 

2) the determination of flushing characteristics between 
basins with multiple connections to inlets 

Consideration of the typical barrier island system indicates 
the systems are comprised of channels, subtidal shallows, tidal flats 
exposed at low tides, and marshes of varying elevation. As the tide 
rises from low to high stages these morphological elements are se- 
quentially flooded. If there is appreciable marsh area the tidal 

N i 

volume increases very rapidly once the marshes become flooded. 

Existing topographic maps and hydrographic maps do not have 
the necessary vertical resolution to obtain accurate estimates of 

the tidal volume as a function of tide elevation. To obtain adequate 

' 1 

coverage using hydrographic and land surveying techniques of the 
tidal flats and marshes would be very expensive and time consuming. 

II The Approach 

The approach used here was to use the water surface itself as 
a contouring machine as the water surface rises over the variable 
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topography with rising tide. This was achieved using sequential 
aerial photo coverage utilising black and white infrared filrn which 
enhances the contrast between flooded and exposed surfaces. In 
order to determine the vertical changes in tidal elevations during 
the sequential "flights "tide gages were installed within the study 
area. The volume storage as a function of tide stage could then 
be computed as the product of incremental increase in the surface 
area flooded and the corresponding incremental rise in tide elevation. 
The details of the procedure are discussed later. 

The experimental test site chosen was that segment of the barrier 
island marsh lagoon complex on the Eastern Shore of Virginia which is 
flooded by the prism gated by Wachapreague Inlet. This study area is 
shown in Figure 25. The rational for this choice was that a substantial 
understanding of the hydraulics of that system had been achieved 
(Byrne, et al., 1973). This background information included direct 
measurements of the prism gated by Wachapreague Inlet and an estimated 
storage relationship based on observation of the tidal stages which 
appeared to initiate flooding of the bays and marshes and a determin- 
ation of the relative area of marshes, bays and channels. The rather 

i * 

poor correlation between measured prisms by integration of instantaneous 
discharge at the inlet and those calculated from the estimated storage 
function led to the decision to attempt the method herein discussed. 

The goal of the study, then, was to obtain a storage curve which 
relates the volume of water within the system serviced by Wachapreague 
Inlet to the tidal stage at Wachapreague Dock tide gage. The tidal - , 
prism for any given tide can then be obtained by subtraction of the 
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low water volume from the high water volume. This capability is 
important to the detailed understanding of inlet response because 
two tides of the same range but differing maximum stage will have 
different prisms and therefore different scouring capacities at the 
inlet channel. 

Ill Methods 

A. Mission Plan. The mission plan was to image the flooded 
surface of the study area at 30 minute intervals as the tide level 
rose from low to high water. During the overflights tidal elevations 
were monitored at eleven sites within the system, and tidal discharges 
were measured at the two interior channels, The Swash and Teagles Ditch, 
which connect the Wachapreague marsh lagoon system with the adjacent 
systems serviced by the adjacent ocean inlets. 

Mission W226, Flight 1, was accomplished on 28 June, 1973, 
utilizing the Wallops Station C54 aircraft equipped with two -T-ll 
aerial mapping cameras using 152mm lenses. The film was typo 2424 
black and white infra-red with a clear anti- vignetting and 89B filter 
combination. The flight was made in hazy weather with broken clouds. 
Visibility was from 5-7 miles. Nominal flight altitude was 9,500 
feet with a wind of 19 knots from 220° . 

Twelve groups of three lines each were flown at approximately 
thirty minute intervals from 1300 to 1900 EDST. The orientation of ' 
the lines made the photography generally suitable for mosaicing; how- 
ever, four of the twelve groups could not be used due to one pr more 
of a combination of excessive tip and tilt, crabbing or lack of 
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sufficient overlap. The complete mission summary is available from 
NASA Wallops Station, Chesapeake Bay Ecological Program Office. 

The flight date was chosen so that the existing conditions 
represented a rising spring tide during daylight hours. A rising 
condition was specifically chosen so that initially drained tidal 
flat faces would be as n dry tT as possible* 

B. Analysis of Data . 

1 ) Ground measurements * The positions of the tide elevation 
measuring stations are shown in Figure 25 and listed m Table 4. 
Stations 5 and 6 are long-term recording tide gages for wharf) tidal 
elevation data planes have been determined. The remainder were 
temporary tide staffs installed on the mission date. Inspection of 
Table 4 indicates a phase lag of up to 0.7 hr. and a range difference 
of up to 0.4 ft. across the system. The most extreme comparison, 
that between Wachapreague Inlet and Wachapreague Dock, is shown in 
Figure 26 referred to long-term mean tide level which is assumed to 
be a level surface. Inspection indicates a water slope of 0.5 - 0.6 
ft. exists between the two positions during the central portion of the 
rising phase'. This fact that the water surfaces slope during the tide 
cycle necessitates the introduction of a correction factor since the 
tide gage at Wachapreague Dock will underestimate the volume, of water 
in the system on the rising tide. Long-term comparisons between these 
two gages show an average range difference of 0.2 ft. and 0.6 hr. 

phase lag* 

Discharge measurements at The Swash (#T1 ir Figure 25) and 
Teagles Ditch (#1 in Figure 25) indicated that tidal water enters 
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TABLE 4 


Station 

1 

2 

3 

4 

5 

6 

7 

8 


9 

. 5.7 

10 

5.6 

11 

5.6 


High Water 
18.5 


11 » 8 hr 

18.4 

11.8 

18.5 

11.3 

18.2 

12.0 

18.7 

11.3 

18.3 

11.4 

18.2 

11.6 

18.5 

11.7 

— 

11.0 

18.3 

12.0 

18.5 


Tide Elevation Data for 28 June, 1973 

Time 

Tide Range (ft.) Low Water 

5.4 

5.5 
5.4 
5.7 
5.3 

5.6 
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WACHA PREAGUE DOCK STAFF 
DEPTH (ft.) 



28 JUNE, 1973 


Figure 26. Tidal elevation curves for Wachapreague Dock and Inlet during 
Mission W 226, Fit. 1. 


the Wachapreague Inlet Basin system during the first part of rising 
tide and exits during the latter part. Measurements on 28 June 1973 
give : 

Swash Teagles Ditch 

Into system 69,000 m 3 Into system 400,000m 3 

Out of system 99,700 m 3 Out of system 548,000m 3 

This will be shown to be a negligible fraction of the total 
volume held in the system. Thus, the effects of this "leakage" are 
ignored hereafter. 

2. Imagery Analysis . Reproductions of the original trans- 
parencies were used in the analysis. Eight of the twelve flight 
groups were assembled into mosaics of the 28.9 square nautical mile 
(86.65 x 10 6 m 2 ) area. 

The major absolute errors in the mosaicing procedure occur in 
aligning some sixteen frames having various edge distortions and in 
determining and cutting conjugate boundaries properly. These errors 
are indeterminate and are presumed random. 

Determinations of the non- water covered areas in each mosaic 
were determined using the J 2 S Digicol located at NASA-Langley . Water 
areas were determined by subtraction from total area- The instrument 
was adjusted so that exposed surface area could be coded as one color 
while the water covered surfaces and the masked areas registered an 
off-scale black. Absolute area calibrations were accomplished using 
masks of known area and plotting these areas versus Digicol- indicated 
percentage area. All mosaics were analyzed on the same day, a con- 
siderable saving in time over hand planimetry. 
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In order to specify film density boundaries associated with 
water boundaries, the water boundaries were visually identified with 
a sharp tipped pointer and the instrument settings were then adjusted 
to shift the boundaries to the proper position. This procedure was 
clear cut at tide levels below the marsh surface. At higher levels, 
marsh plants significantly decreased boundary contrast but sun glint 
assisted in boundary interpretation. For several runs analyzed more 
•than once on different days total water areas varied as much as 5 
percent which revealed a day to day bias in interpretation. These 
problems would likely be reduced if initial film exposure could be 
■adjusted to give high contrast for very shallow depths. The exposures 
during the mission were adjusted to give maximum definition of water 
color differences and not maximum contrast between flooded and non- 
flooded surface. Contrast was heightened in reproduction of the trans- 
parencies. An example of the Imagery is shown in Figure 27a, b for 
a conjugate area at a)low and b) high tide elevations- 

IV Results: Analysis for flooded area and storage volume. 

i 

The times of the flight groups used in the analysis, relative 
to tidal stage are shown in Figure 26. Using flight times and the 
water area determined from the sequential mosaics the relationship 
between . flooded area and stage at the Wachapreague Dock was constructed 
as shown in Figure 28. Observed area changes extended over the range 
in stage from 1.8 ft. to 7.2 ft. at the Dock gage. Since the water 
surface slopes up to as much as 0.6 ft. between the Wachapreague Dock 
and Inlet gages the flooded area measured is not the same as would be 
indicated if the surface was level across the area. 
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Figure 27. Examples of imagery for conjugate areas; a) Flight 
Group 1 and b) Flight Group 12 (time key relative 
to tidal heights is given in Figure 26). 
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Figure 28. Area of flooded survace as a function of tidal stage at 

Wachapreague Dock. Solid line is raw data, dashed line is 
rectified to account for. water slope.. 
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In order to make a first order correction it was assumed that 


at any particular tide stage the surface area flooded is equal to that 
at the stage %AZ (Figure 26) lower. This shifts the curve to lower 
flooded area for a particular stage (dotted line, Figure 28)- The 
segment of the curve between stages 2.5 and 4.3 ft. indicates constant 
flooded area. This represents the near vertical banks of much of the 
major embayment, and that area, 59.5 x 10 s m 2 , represents the cut-off 
between areas of channels, flats, and bays and areas of marshes. This 
and other ancillary information are useful in extending the curve of 
flooded area versus stage to stages lower than the 1.8 ft. observed. 
Visual observations and scattered measurements on the ground indicated 
about 0.5 ft. water depth in the bays at 1130 hrs. on the 28 June. 

Thus, the cut-off for initiation of flooding in the bays was taken 
at 1.3 ft. on the Dock gage. Previous work (Byrne et al., 1973) utilized 
an estimated storage function based upon the areas of channels, marshes 
and bays and the artificial assumption that the -bottom elevations of 
the bays and the marshes were uniform. The measured areas were: 

Channel (major only) 8.087 x 10 e m 2 

Bays 34.8 x 10 6 m 2 

Marsh (including drainage channels) 49.304 x 10 6 m 2 

The cut-off elevations were taken as 2.5 ft. and 5.5 ft. respectively 
for bays and marshes. The curve in Figure 28 indicates that the total 
channel area, including the drainage channel in the marsh, and the 
bays totals to 59.5 x 10 s m 2 . It is reasonable to assume that the 
area of bays is correct at 34.8 x 10 6 m 3 . Thus this total channel 
area must be about 24 x 10 s m 3 . This allows a statement of the component 
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areas within the imaged area as: 

Channels ' 24.0 x 10 s m 2 

Bays 34.8 x 10 6 m 3 

Marsh 37.6 x 10 6 m 3 

Given the .above , the curve in Figure 29 may be extended to the channel 
cut-off; 24 x 10 s m 2 at 1.3 ft. on the Dock gage. 

The cumulative volume of water stored in the marsh- bay- channel 
system was computed from the curve in Figure 28 using incremental 
volumes for every 0.5 ft. stage increase at the Dock gage. The volumes 
are referred to 0 elevation of that gage which insures inclusion of 
all observed tides within the last three years. The assumption was 
made in the calculations that the flooding for each, incremental in- 
crease in stage occurred over a uniformly sloping bottom. The cumula- 
tive volume versus gage is shown in Figure 29. The figure also shows 
the originally estimated volume using channel, bays and marshes as 
components with discreet cut-off elevations. The curve has been extended 
from the 7.2 ft. to the 10 ft. stage since visual observation indi- 
cated the entire system was flooded at stages beyond 7.2 ft. The 
shift to the left of the storage curve based on the imagery relative 
to that originally estimated is due to adoption of the 1.3 ft. stage 
as the cut-off between channel flooding and the initiation of bay 
flooding . 

The storage relationships shown in Figure 29 are compared with 
tidal volumes measured at Wachapreague Inlet in Table 5'. The observed 
values were determined by integration of instantaneous tidal discharge 
measured by current meter arrays at the inlet throat. Although 
accuracy of the directly measured prisms is not known, the measurements 
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Figure 29. Storage function for Wachapreague Inlet basin as determined 
by infra-red imagery and as originally estimated. 
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TABLE 5 


Comparison of measured tidal volumes at Wachapreague Inlet 
with volumes indicated by imagery and estimated storage. 


TIDAL PRISM (m 3 ) % Dif f (3+1) 


Tide 

Cycle 

Number 


Date 

Tide stage at 
high and low 
m (ft.) 

Po - 

measured 
at inlet 

Pi 

from 

imagery 

Ps 

from 

estimated 

i 

?! 

P S 

1 

14 

July 70 

2. 01 (6*6)-. 79 (2. 6)' 

64.9 x 10 6 

78 x 10 s 

67.9 

x 10 s 

+20.2 

+ 4.7 

2 

14 

July 70 

2.01(6.6)- .85(2.8) 

62.2 

74.5 

65.3 


+19.8 

+ 5.0 

3 

25 

Aug. 71 

1.86(6*1)“. 79(2. 6) 

57.4 

66 

66 

53.9 


+15.0 

-'6.2 

4 

13 

Sept. 72 

1.95(6.4)-. 98(3. 2) 

63.2 

63 

54.4 


~ 0.3 

-13.9 

5 

13 

Sept . 72 

1.73(5. 7)-. 98(3. 2) 

47.6 

47.5 

34. 7 


- 0.2 

-27.0 

6 

14 

Sept . 72 

1.73(5. 7)-. 94(3.1) 

45.5 

49 

36 


+ 7.7 

-20.8 

7 

14 

Sept . 72 

1.86(6.1)-. 94(3.1) 

61.0 

58 

47.3 

i 

- 4.9 

; 

-22.5 






represent the best available estimates. Comparison of cycles 3 and 7 
in Table 5 indicates an inconsistency in the measured values; both 
cycles had the same maximum stage but the minimum stage of cycle 3 
was 0.15m less than cycle 7. In actuality cycle 3 must have had a 
larger prism than cycle 7 but the measurements indicate the converse. 

A similar inconsistency exists between cycles 5 and 6 although to a 
lesser degree. The data for measured prisms indicated the accuracy 
is no better than ±10%. 

Inspection of the percent differences between measured prism 
and that deduced from the imagery and the original estimated storage 
indicates that volumes deduced from the imagery are in appreciably 
better agreement with the measured than are those from the estimated 
storage. The average absolute differences are 9.7% and 14.3% re- 
spectively.' As previously mentioned the "leakage” from the Wacha- 
preague storage system via the interior channels, The Swash and 
Teagles Ditch, does not appear to significantly influence the error 
in estimating the storage from the imagery; the prism for the tide 
cycle used to construct the storage curve was 106.5 x 10 6 m 3 , while 
data listed in earlier paragraphs for the volume passage in The Swash 
and Teagles Ditch combined amount to less than 1% of that figure. 

Thus, it appears justified to disregard this contribution to the storage. 

V Discussion 

The present study demonstrates the feasibility of utilizing 
sequential black and white infra-red photography to determine the 
storage characteristics of complex basins wherein the flooded area 
changes irregularly with varying stages of rising and falling tides. 
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This procedure resolves the problems encountered with determination of 
tidal prism based on estimates of basin area at high and low water 
due to the ambiguity of these determinations in shallow irregular 
basins (0 T Brien, 1969; 0 T Brien and Clark, 1973). The method offers 
great potential value 'in tidal inlet studies as it offers the op- 
portunity to specify the expected average discharge for any particular 
tidal condition. Furthermore, if the inlet channel cross-sectional 
area is known, the average ebb and flood currents for any given tidal 
cycle may be derived. This information may be then summed over a 
longer period (months, seasons, years) for which reference tide records 
exist to estimate whether the inlet currents advect sand into or out 
of the inlet (Byrne et al. , 1973). 

The methodology also offers very interesting possibilities as 
an alternate approach to constructing depth contour maps. The primary 
mission of the National Ocean Survey *s hydorgraphic mapping activities 
is to provide highly reliable maps for navigation purposes. Typically, 
in areas where such studies are done the channels are carefully delineated 
but less emphasis is given to the shallow bays and no attention is 
given to the variable marsh elevation. State of the art numerical ' 
models for flushing such as the Jamaica Bay Study (Leendertse et al. , 

1973) require moving boundary conditions which necessitates fairly 
high resolution depth maps. Such maps could be constructed using the 
described technique wherein the instantaneous waterline acts as a 
contour. This approach would be quite successful in those areas 
where there is not appreciable tide range variation within the basin 
system; such as the Wachapreague System. 
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Finally, it should be noted that the technique offers a pos- 
sible method to study tidal range variations in those basins with large 
range and phase differences within the system. Existing methods in- 
volve the installation of a large number of tide gages and a data base 
over a period of one month or longer. The tidal datum. plane is then 
established and mean range and phase differences determined. In those 
areas where such detailed surveys cannot be conducted the method dis- 
cussed herein would supply very useful data. It is also likely that 
application of the method would be useful in planning tide. gage lo- 
cations for a detailed survey. In such applications the stage variations 
as a function of time would be measured using stereophotogrammetric 
techniques. 
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11.7 

10.4 

9.2 

9.2 

7.1 

5.3 

12.8 

4 44* 

* * * * 

4444 

10.5 

26 

11.1 

8.6 

5.5 

4.4 

3.2 

3.3 

6.1 

6.6 

9.5 

4 4 4 4 

aa a* 

4444 

6.5 

9 

14.0 

13.2 

11.7 

11.1 

10.4 

9.8 

8.1- 

7.8 

12.1 

4 4 4 * 


444 4 

10.9 


FLIGHT 2 RUN 3 


POINT 

FR 1 

FR 2 

FR 3 

11 

3.7 

2.9 

1-5 

Ml 

6.7 

6.7 

6.2 

33 

2.1 

2.1 

3.5 

25 

3,0 

2.0 

1.6 

16 

B .6 

8.1 

5.3 

19 

6.6 

3.3 

2.0 

M2 

3.8 

3.7 

6.7 

15 

6.9 

6.9 

6.9 

18 

3.5 

2 .9 

2.1 

26 

4**4 

13.6 

8.6 

10 

4 4*4 

*4*4 

9.8 

22 

* 4 * 4 

*44* 

9.8 

26 

10.7 

10.1 

7.3 

20 

9.8 

8.6 

5. 1 

9 

10.6 

9.0 

8.6 


FLIGHT 

.SUMMARY, 

FLIGHT 

POINT 

AVG 

S ICM A 

11 

6.63 

2.12 

Ml 

6 . 30 

0.78 

33 

3.67 

1.75 

25 

6.66 

2.0.’ 

16 

7.77 

2.96 

19 

3.60 

2. CO 

M2 

5.32 

2.89 

15 

6.31 

2.86 

18 

3.61 

1.83 

26 

11 .95 

6.33 

10 

7.69 

3.66 

22 

8.15 

3.23 

26 

10,63 

3.39 

28 

6. 98 

2.51 

9 

11.15 

2.8 3 


FR 6 

REFLECTANCES 
FR 5 FR 6 FR 7 FR 8 

FR 9 

FRIO 

FRil 

FR12 

AVG 

2.0 

2.1 

2.9 

3.7 

6.1 

A*** 

A * A A 

AAA* 

AAA* 

3.1 

6.5 

3. 7 

6.2 

5.1 

6. 1 

AA* A 

AAA * 

A A A A 

A A A A 

4.7 

3.8 

6.1 

44*4 

* 44 * 

4 444 

A * * * 

* * A A 

AAA* 

A * A A 

3.5 

2.2 

2.9 

2.2 

3.8 

44*4 

*4 A A 

A * * A 

A AAA 

AAA* 

2.5 

5.7 

3.8 

3.2 

6.2 

6.9 

As*** 

A AAA 

A A A A 

AAA A 

5.4 

0,9 

0.9 

1.3 

2.3 

3.8 

*** * 

A A A A 

A A A A 

A AAA 

2.4 

6.9 

fc. 1 

7.6 

44*4 

44* * 

A*** 

A A A A 

A A A A 

A A A * 

5.1 

5.5 

8.6 

13.6 

4444 

*444 

A A A * 

A A A A 

AAA A 

A* A A 

7.0 

1.7 

2.C 

2.6 

6.0 

6. 1 

*4* A 

A AAA 

A A A A 

A * A* 

3.1 

8.1 

6.7 

4.7 

6.4 

8.9 

12.1 

A A A A 

A A A A 

A A A A 

B.4 

7.3 

3.8 

6.5 

2.9 

5.1 

S.7 

A A A A 

A A A A 

AAA* 

5.6 

8.9 

5.1 

3.0 

3.5 

6.5 

5 4 3 

A A A A 

A A A A 

AAA * 

5.9 

7.6 

5.9 

7.3 

7.1 

7.1 

9 * 2 

A AAA 

A A A A 

A A A A 

e.o 

6.1 

5,7 

5.7 

7.3 

9.5 

12-5 

AAA A 

AAA A 

A A A A 

7.8 

10.7 

7.1 

7.L 

7.1 

8.9 


* A A A 

A A A A 

* A A A 

8.7 


2 


gglODUCIBILITY OS' THE 
CSMNAL PAGE IS POOR 


» ( 
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SIGMA 

L.8 

0.7 

2.2 

1.4 

1.9 

1.8 

4.6 
1.1 
1.8 

4.2 

3.7 

2.7 

2.1 

2.3 
2.3 


SIGMA 

1.4 
0.7 
1 .0 
1.6 
2.6 
2.0 
0.3 
1.6 
1.2 
2.6 

2.3 

2.5 

3.3 
2.7 
2.1 


SIGMA 

1.4 

0.7 

1.6 

0.7 

1.8 

1.3 

1.4 

3.4 

1.4 

3.2 

2.3 

2.4 

1.5 

2.4 

1.5 
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DATE 

STARTING 

SLACK 

TRANSECT 





MO DA 

YR 

TIME 

WATER 

NUMBER 





4 17 

73 

11 2 

LOW 

4 






WATER 

TEMPERATURE 

AIR 

TEMP. 

DEW PT. 

SALINITY 


TIME 

0.5FT 

3FT 

6FT 

3FT 

6FT 

TEMP. 

0/00 

HR 

MIN 

SEC 

D.F 

D.F 

D.F 

D.F 

D.F 

D.F 


11. 

2 

36 

56 • 54 

56.38 

56.46 

67.25 

66.91 

66.25 

2.00 

11 

3 

1_2 

56.71. 

56.54 

56.62 

67.50 

66.75 

51.90 

2.03 

11 

3 

48 

56.87 

5 6.71 

56.79 

67.25 

6 6.75 

52.83 

2.11 

11 

4 

24 , 

57.12 

56.96 

57.04 

67.16 

66.42 

53.24 

1.98 

11 

5' 

0 

57.46 

57.21 

57.29 

67.25 

66.42 

53.37 

2.07 

11 

5 

36 

57.62 

57.46 

57.46 

67.25 

66. 50 

53.10 

2.03 

11 

6 

12 

57.79 

57.70 

57.70 

67.41 

66.25 

52-70 

2.06 

11 

6 

48 

57.95 

57.87 

57.79 

67.16 

66.58 

53.77 

1.90 

11 

7 

24 

58.04 

57.87 

57.79 

67.33 

66.83 

52.30 

2.10 

11 

8 

0 

58.04 

57.87 

57.62 

67.08 

66.83 

53.10 

2.05 

11 

8 

36 

58.04 

57.87 

57.46 

67.41 

67.08 

52.43 

2.06 

11 

9 

12 

57.95 

5 7.79 

57.70 

67.75 

67. 33 

51.64 

2.05 

11 

9 

48 

57.79 

57.79 

57.95 

67. 50 

67.25 

52.70 

2.13 


73 






DATE 

STARTING 

SLACK 

TRANSECT 





MO DA 

YR 

TIME 

WATER 

NUMBER 





4 17 

73 

1110 

LOW 

5 






WATER 

TEMPERATURE 

AIR 

TEMP. 

DEW PT. 

SALINITY 


TIME 

0.5FT 

3FT 

6FT 

3FT 

6FT 

TEMP. 

0/00 

HR 

MIN 

SEC 

D.F 

D.F 

D.F 

D.F 

D.F 

D.F 


11 

10 

32 

58.78 

58.54 

58.62 

67.58 

66.91 

53.64 

2.04 

11 

11 

4 

58. 04 

57.87 

57.87 

67.50 

66.50 

53.37 

2.16 

11 

11 

36 

57.46 

57.46 

57.46 

6 7 . 58 

66.50 

52.57 

2.05 

11 

12 

a 

57.54 

57.46 

57.12 

67.66 

66.58 

52.57 

2.07 

11 

12 

40 

57.70 

57.54 

57.21 

67.75 

66.00 

53.64 

1.95 

11 

13 

12 

57.87 

57.70 

57.29 

67.83 

66* 50 

52.43 

1.94 

11 

13 

44 

57.79 

57.54 

57.46 

68.00 

66.58 

52.83 

1.96 

11 

14 

16 

57.79 

57.54 

57.54 

68.25 

67.08 

52.17 

2.06 

11 

14 

48 

57.79 

57.62 

57.62 

68.08 

66.66 

52.83 

1.99 

11 

15 

20 

57.62 

5 7.79 

57.54 

68.25 

67.00 

52.17 

1.96 

11 

15 

52 

57.62 

57.46 

57.46 

68.16 

66. 83 

52.57 

1.95 

11 

16 

24 

57.54 

57.46 

57.37 

68.33 

66.75 

52.70 

2.01 

11 

16 

56 

57.46 

57.37 

57.37 

68.25 

66.66 

52.70 

1.98 

11 

17 

28 

57.37 

57.29 

57.37 

68.33 

66.83 

52.57 

1.96 

11 

18 

0 

57.12 

57.04 

57.12 

68.41 

66.91 

52.04 

2.05 

11 

18 

32 

56.87 

56.79 

56.96 

68.33 

67.08 

53.24 

1.94 

11 

19 

4 

56.87 

56.71 

56.87 

68.75 

67.50 

51.90 

1.99 

11 

19 

36 

57.12 

57.04 

57.04 

68.58 

68.16 

50.85 

2.04 

11 

20 

8 

57.54 

57.46 

57.37 

68.58 

67.33 

53.10 

1.97 

11 

20 

40 

5 7.62 

57.54 

57.54 

68. 58 

66. 83 

53.24 

1.95 

11 

21 

12 

57.46 

57.21 

57.21 

68.91 

67.08 

52.43 

2.06 

11 

21 

44 

56.71 

56.62 

56.62 

68.91 

67.58 

51.77 

2.12 

11 

22 

16 

56.71 

56.54 

56.46 

69.00 

68.41 

51.77 

1.93 

11 

22 

48 

56.46 

56.29 

56.21 

70.08 

68.33 

51.24 

2.03 

11 

23 

20 

56.29 

56.13 

56.04 

69.25 

67.91 

51.64 

1.94 

11 

23 

52 

56.13 

55.96 

55.96 

69. 16 

67.91 

51.64 

2.03 

11 

24 

24 

56.21 

56.04 

56.04 

69.25 

67. 33 

52.70 

1.95 

11 

24 

56 

56.04 

55.96 

55.88 

69.08 

68.00 

51.24 

1.90 

11 

25 

28 

55.96 

55.79 

55.54 

69.08 

67.83 

51.51 

1.95 

11 

26 

0 

55.63 

55.13 

55.04 

69.42 

68.66 

50.07 

1.93 

11 

26 

32 

55.54 

55.21 

55.21 

69.25 

67. 75 

50.98 

1.94 

11 

27 

4 

55.63 

55.13 

55.04 

69.00 

66.83 

51.24 

1 • 90 

11 

27 

36 

55.63 

55.38 

55.21 

69.00 

67.50 

50.85 

i .98 

11 

28 

8 

55.54 

55.21 

55.21 

69.42 

68. 00 

50.46 

1.81 

11 

28 

40 

55.21 

55.04 

55.13 

69.00 

67.33 

51.38 

2.03 

11 

29 

12 

54.96 

54.79 

54.88 

69.25 

67.00 

50.72 

1.93 

11 

29 

44 

54.96 

54.96 

54.96 

69.00 

66.91 

51.90 

1.93 

11 

30 

16 

55.21 

55.04 

55.13 

69.42 

68.25 

50.46 

1.92 

11 

30 

46 

55.13 

55.04 

55.13 

69.42 

68.58 

50.59 

1.90 
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DATE 

STARTING 

SLACK 

TRANSECT 





MO DA 

VR 

TIME 

WA TER 

NUMBER 





4 17 

73 

1132 

LOW 

6 






WATER 

TEMPERATURE 

AIR 

TEMP. 

DEW PT. 

SALINITY 


TIME 

0.5FT 

3FT 

6FT 

3FT 

6FT 

TEMP. 

0/00 

HR 

MIN 

SEC 

D.F 

D.F 

D.F 

D.F 

D.F 

D.F 


11 

32 

34 

55.13 

55.04 

55.04 

69.58 

68.83 

51.64 

2.01 

11 

3-3 

-8- - 

55.21 

55-. 0-4- 

5-5-. 1 3 

69.58 

69.33 

51.77 

1.89 

11 

33 

42 

55.13 

55.13 

55.13 

69.50 

68.83 

52.17 

1.94 

11 

34 

16 

54.88 

54.96 

54.79 

69.33 

68.50 

52.57 

1.91 

11 

34 

50 

55.21 

54.96 

55.04 

69.16 

68.33 

52.43 

1.92 

11 

35 

24 

55.29 

55.04 

54.96 

68.91 

68.08 

52.97 

1.89 

11 

35 

58 

- 54 * 96 

54.88 

54.79 

69.33 

68.33 

52.17- 

1.92 

11 

36 

32 

55.21 

55.21 

54.96 

69.08 

68.00 

52.70 

1.94 

11 

37 

6 

55.96 

55.71 

55.29 

69.16 

68.33 

52.97 

1.93 

11 

37 

40 

56.04 

55.71 

55.38 

69.08 

68.08 

52.83 

1.94 

11 

38 

14 

56.13 

55.29 

55.21 

69. 16 

68.08 

52.70 

2.03 

11 

38 

48 

56.04 

55.54 

55.38 

69.08 

68.16 

52.97 

1.93 

11 

39 

22 

56.21 

55.71 

55.46 

69.08 

68.08 

52.97 

1.87 

11 

39 

56 

56.29 

55.71 

55.54 

69.33 

68.41 

52.57 

1.85 

11 

40 

30 

56.46 

56.13 

55.71 

69.16 

68.41 

52.17 

1.92 

11 

41 

4 

56.71 

56.21 

55.96 

68.91 

68.41 

52.70 

2.02 

11 

41 

38 

57.04 

56.54 

56.04 

69.33 

68.58 

52.04 

2.01 

11 

42 

12 

57.46 

57.04 

56.13 

69.50 

68.08 

52.43 

1.89 

11 

42 

46 

57.79 

56.79 

56.21 

69.58 

68.08 

52.57 

1.89 

11 

43 

20 

57.70 

57.04 

56.21 

69.50 

68.08 

52.43 

1.89 

11 

43 

54 

57.70 

56.87 

56.38 

69.58 

68.50 

53.37 

1.94 

11 

44 

28 

57.87 

57.12 

56.71 

69. 33 

68.50 

53.24 

1.93 

11 

45 

2 

57.87 

57.37 

56.79 

69.25 

68.08 

53.37 

1.88 

11 

45 

36 

57.79 

57.37 

56.62 

69.08 

67.91 

53. 64 

1.94 

11 

46 

10 

57.54 

57.21 

56.54 

69.08 

67.91 

53.50 

1.93 

11 

46 

44 

57.37 

57.04 

56.29 

69.00 

68. 00 

53.24 

1.93 

11 

47 

18 

57.37 

57.04 

56.04 

69.00 

67. 91 

52.97 

1.95 

11 

47 

52 

57.29 

57.12 

56.54 

68.83 

67.25 

53.37 

1.92 


75 



DATE STARTING SLACK TRANSECT 





MO DA 

YR 

TIME 

WATER 

NUMBER 





4 17 

73 

1149 

LOW 

7 






WATER 

TEMPERATURE 

AIR 

TEMP. 

DEW PT. 

SALINITY 


T IME 

0.5FT 

3FT 

6FT 

3FT 

6FT 

TEMP. 

0/00 

HR 

MIN 

SEC 

D • E 

D.F 

D.F 

D.F 

D.F 

D.F 

11 

49 

30 

56.96 

56.87 

56.29 

69.33 

68.41 

51.51 

1.98 

11 

50 

0 

57.04 

56.96 

55.88 

69.67 

68.83 

51.11 

2 a 03 

11 

50 

30 

57.04 

5 6. 46 

55.71 

69.25 

68.16 

52.04 

1.88 

11 

51 

0 

56.96 

56.38 

55.63 

69.25 

67.41 

52.97 

1.88 

11 

51 

30 

56.79 

56.62 

55.46 

69.25 

67.75 

52.17 

1.95 

11 

52 

0 

56.71 

56.46 

55.21 

69.16 

68.08 

52.04 

1.92 

11 

52 

30 

56.29 

56.13 

55.04 

69.33 

67. 91 

51.90 

1.89 

11 

53 

0 

55.96 

55.79 

54.79 

69.42 

67.83 

52.57 

1.95 

11 

53 

30 

56.46 

56.04 

54.54 

69,00 

67.75 

52.83 

1.93 

11 

54 

0 

56.46 

55.38 

54.46 

69.25 

67.58 

51.51 

1.92 

11 

54 

30 

54.54 

54.46 

54.46 

69.42 

68.08 

51.38 

1.91 

11 

55 

0 

54.54 

54.46 

54. 46 

69.67 

68. 66 

50.98 

1.92 

11 

55 

30 

54.63 

54.71 

54.71 " 

69.92 

68.25 

51.51 

1.88 

11 

56 

0 

57.04 

56.71 

55.04 

69.67 

68.66 

52.17 

1.91 

11 

56 

30 

57.70 

57.12 

55.21 

69.58 

68.16 

53.10 

1.89 

11 

57 

0 

58.12 

57.54 

55.46 

69.67 

68.91 

51.90 

1.89 

11 

57 

30 

58.04 

57.70 

56.04 

70.08 

69. 00 

51.24 

1.90 

11 

58 

0 

58.20 

57.70 

56,79 

69.92 

69. 16 

51.51 

1.97 

11 

58 

30 

58.29 

57.70 

56.87 

70.00 

69. 00 

51.77 

1.99 

11 

59 

0 

58.54 

58.04 

56.87 

70.17 

68.83 

51.24 

1.91 

11 

59 

30 

58.78 

58.37 

56.71 

70.33 

69.50 

51.11 

1.92 

12 

0 

0 

58.87 

58.37 

56.62 

70.17 

69.25 

51.38 

1.65 

12 

0 

30 

58.78 

58.12 

56.29 

69.83 

68. 33 

51.90 

1.65 

12 

1 

0 

58.37 

57. 79 

56.04 

69.63 

68.58 

51 . 64 

1.65 

12 

1 

30 

57.62 

57.12 

55.79 

69.83 

68.75 

51.24 

1.65 

12 

2 

0 

56.79 

56.71 

55.29 

69,58 

68.41 

52.04 

1.65 

12 

2 

30 

55.46 

55.21 

55.13 

69.75 

68.66 

51.38 

1.65 
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DATE STARTING 
MO DA YR TIME 

A 17 73 12 3 





WATER 

TEMPERATURE 


TIME 

0.5FT 

3FT 

6FT 

HR 

MIN 

SEC 

D.F 

D.F 

D.F 

12 

3 

20 

56.13 

56.13 

55.21 

12 

3 

40 

57.54 

-5-7— 0-4- 

55.54 

12 

A 

0 

58.20 

57.54 

55.63 

12 

A 

20 

58.78 

57.87 

55.88 

12 

4 

40 

58.87 

57.79 

56.04 

12 

5 

0 

58.95 

58.20 

56.38 

12 

5 

20 

50.95 

58.37 

56.29 

12 

5 

40 

58.62 

58.29 

56.54 

12 

6 

0 

58.29 

57.46 

56.62 

12 

6 

20 

58.37 

57.95 

56.13 

12 

6 

40 

58.20 

57.87 

56.62 

12 

7 

0 

58.20 

57.79 

56.29 

12 

7 

20 

57.87 

57.54 

56.29 

12 

7 

40 

55.21 

54.79 

54.71 

12 

8 

0 

54.88 

54.71 

54.71 

12 

8 

20 

54.88 

54.71 

54.71 

12 

8 

40 

54.88 

54.79 

55.21 

12 

9 

0 

56.46 

55.88 

55.13 

12 

9 

20 

56.13 

55.38 

55.04 

12 

9 

40 

55.79 

55.30 

55.21 

12 

10 

0 

55.79 

55.71 

55.29 

12 

10 

20 

55.88 

55.29 

55.13 

12 

10 

40 

55.79 

55. 71 

55.63 

12 

11 

0 

55.71 

55.63 

55.71 


SLACK TRANSECT 
WATER NUMBER 
LOW 8 


AIR 

TEMP. 

DEW PT. 

SALINITY 

3FT 

6FT 

TEMP. 

0/00 

D.F 

D.F 

D.F 


70.25 

69.67 

50.07 

1.65 

70.42 

69c 67 

- 51.3.8 

1.65 

70.33 

69.58 

51.64 

1.65 

70.17 

69.25 

52.70 

1.65 

70.50 

69.75 

51.24 

1.65 

68.33 

69.75 

52.43 

1.65 

70.59 

70.25 

51.51 

1.65 

70.50 

70.25 

51.90 

1.65 

70.75 

70.50 

50. 72 

1.65 

70.84 

70.25 

51.24 

1.65 

70.92 

70.25 

51.51 

1.65 

71.00 

70.42 

51.38 

1.65 

70.92 

70.25 

52.04 

1.65 

70.92 

70.25 

52.04 

1.65 

71.00 

70.25 

52.30 

1.65 

71.00 

70.17 

51.90 

1.65 

70.92 

70.08 

51.90 

1.65 

71.17 

70.00 

51.90 

1.65 

70.92 

70.08 

51.64 

1.65 

71.26 

70. 08 

51.51 

1.65 

71.26 

70.08 

51.51 

1.65 

71.17 

70.17 

51.11 

1.65 

71.34 

70.25 

51.11 

1.65 

71.09 

70.08 

52.04 

1.65 
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